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Abstract. The effects of saving and spending patterns on holding time distribution of money are investi- 
gated based on the ideal gas-like models. We show the steady-state distribution obeys an exponential law 
when the saving factor is set uniformly, and a power law when the saving factor is set diversely. The power 
distribution can also be obtained by proposing a new model where the preferential spending behavior is 
considered. The association of the distribution with the probability of money to be exchanged has also 
been discussed. 
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1 Introduction studying of the distribution functions, is essential. Some 

^ ■ 

. pioneering work along this line has been reviewed in a 

$H ' During the last several years, methods and techniques of 

popular article [T5| . 

statistical physics have been successfully applied to eco- 



As well known, the exploration of the distribution of 
money can be traced back at least a century to the work 
of the Italian social economist Vilfredo Pareto, who stud- 
ied the distribution of income among people in different 
western countries and found an inverse power law |14| . 

Recently this topic has been taken up with the emergence 

Send offprint requests to: Yougui Wang 

E-mail: ygwang@bnu.edu.cn of econophysicists among whom some believe that there 



nomical and financial problems 1,2,3 . Recently, Some 
econophysists have been paying attention to the statistical 
mechanics of money, theoretically or empirically | 4l5lfil7l 
8,9, ITTlllllli2| . They believe that a thorough understand- 
ing of the statistical mechanics of the money, especially 
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might be some physical and mathematical rules govern- intervals for money to be held. This kind of time interval 
ing the distribution of income or wealth in the world and was ever called by Wicksell the "average period of idle- 
attempt to discover them. A series of models have been ness" or "interval of rest" of money |15| . In our previous 
developed for the equilibrium money distribution based on work , we called it "holding time" of money and found 
the analogy between market economics and kinetic theory that after the economy has achieved an equilibrium state, 
of gases 4,5 6,812 . Identifying exchange between any there is not only a distribution of money among agents, 
two agents in a closed economy where the money is con- but also a steady distribution over the holding time. We 
served with the two-body elastic collision in an ideal gas, also found that monetary velocity, an important macroe- 
these models show no matter how uniformly and forcefully conomic variable, which is associated with Irving Fisher 
one distributes money among agents initially, the succes- J7| , could be expressed as the expectation of the recipro- 
sive tradings eventually lead to a steady distribution of cal of holding time, 
money. And the shape of money distribution is determined 
by the trading rule for choosing an amount of money to 
transfer. Allowing agents to hold back some of their money 
when they are chosen to trade, B.K. Chakraborti et al. 

In a basic ideal gas-like model, the distribution of money 

introduced the saving behavior into the model by adding 

over the holding time follows an exponential law, where 

a saving factor s in the trading rule UJ. The simulation 

saving behavior is not taken into account. The purpose 

results clearly indicate a robust Gibbs-like distribution 

of this paper is to study how the introduction of saving 

where the density of agents with money m decreases ex- 
behavior affects such kind of distribution. In next section, 

ponentially with m for s — 0, which is identical to the 

we make a brief review of the basic ideal gas-like model 

result of A. Dragulescu and V. M. Yakovenko's random 

by which our work can be erected and of the measure- 
two-agent exchanges model : 5 J. The distribution of money 

ment of the distribution of holding time. In sections and 

changes to follow asymmetric Gibbs-like law when the 

01 we show that the uniform saving factor gives exponen- 

fixed and uniform saving factor is set to be nonzero, while 

tial distribution, while the diverse saving factor induces a 

a 'critical' Pareto distribution of money is found when 

change to power distribution. Then we introduce prefer- 

saving factor is set diversely among agents [E1EIIH]- 

ential spending behavior into the model in section [S] and 
In practice, money is held as a store of value, what again obtain power distribution. Comparing these results, 
is more, it plays an essential role for being a medium of we can conclude that the formation of holding time distri- 
exchange. Money is transferred consecutively from hand bution is associated with the character of the probability 
to hand in the exchange process, in which there exist time of money to be exchanged. 
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2 An Ideal Gas-Like Market Model and 
Holding Time Distribution 

We begin with the basic ideal gas-like model which was 
introduced firstly by A. Dragulescu and V.M. Yakovenko 
A close economy is considered in the model where 
the amount of money M is conserved and the number 
of agents N is fixed. The money is possessed by agents 
individually and agents can exchange money with each 
other. Since the scale and initial distribution of money 
have no effect on the final results, most of our simulations 
were carried out with N = 250 and M = 25000 and the 
amount of money held by each agent was set to be M/N 
at the beginning. The trade in the economy is modelled to 
take place round by round. In each round, two agents, i 
and j for example, are chosen randomly to get engaged in 
a trade among which agent i is "receiver" and the other 
one j is "payer" . The amount of money that changes hand 
Am is determined by trading rule which ensures that the 
amount of any agent's money is non-negative and the total 
money is conserved. A trading rule commonly used can be 
expressed as Am = s(rm + rrij)/2, where e is a random 
number from zero to unity. As for which units of money 
are chosen to be transferred, all in the payer's hand is 
equally probable. 

In the ideal gas-like model, money is held by agents 
and transferred frequently. In this process, if an agent re- 
ceives money from other agents, he will hold it in hand 
till paying it out to some other agents. The time interval 
between the receiving and paying out is named as holding 
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time |lfij . The holding times of a certain unit of money at 
different moments or those of different units of money at a 
given moment are not the same. We introduce the proba- 
bility distribution function of holding time Ph(r), which is 
defined so that the amount of money whose holding time 
lies between r and r + dr is equal to MPh{r)dr. So, we 
can get the normalization condition and the expression of 
the expectation of holding time as follows: 

POO 

/ P h {r)dT = l (1) 
Jo 

and 

/>oo 

T= / rP h (r)dT. (2) 
Jo 

In the simulations, suppose it is at round to that we 
start to reord, and so, holding time is recorded as the 
difference between the moments when the money takes 
part in trade after to for the first two times. The recording 
mode is illustrated in Figure 1(a). Please note this mode 
is different from what we have adopted in Reference |16|. 
which is shown in Figure 1 (b) . The measurement results of 
the two modes seem quite different, however, they reflect 
the same process in different ways. We adopt herein the 
mode (a) solely to facilitate the exposition. The typical 
distribution of holding time is shown in Figure 2. It can 
been seen from the inset of Figure 2 that the distribution 
of holding time follows an exponential law: 

iMr) = ie-*. (3) 

This result indicates that the transferring process of money 
is a Poisson process with intensity of ^. 

To get the distribution of holding time without system- 
atic factor disturbing, we performed the simulations about 
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100 times with different random seeds and data were not and the number of agents is fixed. Any agent's money 
collected until the probability distribution of money got is non-negative or no debt is allowed. The agents are in- 
stationary. And for convenience we stopped data collecting distinctive at the beginning of simulations: same initial 
after majority of money(> 99.9%) had been recorded. In amount of money and same saving factor s. In each round, 
all the following simulations, holding times are measured an arbitrary pair of agents are chosen to make exchange 
in this way after the distributions of money get stationary with each other. For example, at t-th round, agent i and j 
of course. take part in trading, so that at t + 1-th round their money 

rrii(t) and mj(t) change to 

3 Model with Uniform Saving Factor , , 

m;(f + 1) = m t {t) + Am; m,j(t + 1) = m 3 -(f) - Am; (4) 

In reality, saving behavior is a natural action pattern for 

where 

any economic agent. In order to insure future consump- 
tion, people always keep a part of their money as saving. ^ m = — s )i( £ ~ + £m j w]i (5) 
The ratio of the saving to total amount of money held by 

stitution, it is obvious that the trading rule satisfies the 



and e is a random fraction. After a straight-forward sub 

an agent is called "marginal propensity to save" by B. K 



Chakrabarti's group. The term of marginal propensity to . . ,. . .... . . 

conservation and non- negativity condition, and each agent 

save has totally different meaning in economics, which is . r , . . . . 

saves fraction s of his money before trade. 

defined as the partial derivative of saving function with re- _,. . . . . . . _,. „ „ 

Ihe simulation results are shown m figure 6, for some 

spect to income To avoid confusion, we rename it as . „ T , , , .. . . .. , ., 

^* values of s. It can be seen that the probability distribu- 

"saving factor". Referring to the saving factor, two cases . (.•■•,,. • r n n i 

tions of holding time for all saving factors decay cxponcn- 

have been considered by Chakrabarti's group, one is that . . , . . . . , . 

tially. And the lower the saving factor is, the steeper the 

all agents have a uniform saving factor, the other one is .... . .... , . , . , r 

distribution curve. I hese results indicate this kind of sav- 

that saving factors are randomly distributed among agents . , . . , . _ . r , . 

mg behavior does not change the Poisson nature of the 

IBHSI. As mentioned above, they found the equilibrium dis- , . . . ., . , 

exchanging process, but its intensity. 

tributions of money among agents had remarkable differ- 
ent characters under such two assumptions. Along this 

4 Model with Diverse Saving Factor 

line, in this section and the next one we shall examine 

the impacts of the saving behavior on the distribution of In realistic economy, how much an agent saves depends 
holding time for the two cases respectively. on the economic situations he or she faces, and the saving 

All the assumptions of the above ideal gas-like model factor of course varies from agent to agent due to their 
do work in this model. The amount of money is conserved different conditions. To get closer to reality, this model 
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inherits all the assumptions and evolution mechanism of non-linear associations. We also found the correlation co- 

the previous model except that of uniform saving factor, efficient falls to and keeps at about 0.32 after 500000th 

Each agent's saving factor is initialized at the beginning round. Although the value of the correlation coefficient 

of simulations which distributes randomly and uniformly is not high enough, it still implies that the agents with 

within an interval to 1, and is fixed in the simulations, higher saving factors hold more money. 

Correspondingly, the trading rule Equation © changes to Secondly, we computed the average value of saving 



Am = (1 - s t )(e - l)m t {t) + (1 - s^em^t)- (6) 



where Si, Sj are the saving factors of agent i and j respec- 
tively. 

To our surprise, once the diverse saving factor is intro- 
duced into the model, as shown in Figure 4, the holding 
time distribution changes to obey a power law instead of 
an exponential law. This result indicates that the transfer- 
ring process of money in this model is not a Poisson pro- 
cess any more. The Poisson nature of the process is broken 
due to the loss of homogeneity of the money transferring. 
In the previous model, for any saving factor, the proba- 
bility of each unit of money to participate in exchanges 
at any round is equal because the saving factor is set to 
be uniform for all agents. On the contrary, in this model 
the transferring probability of money is not the same any 
more due to the diversity of the saving factors. This con- 
clusion was verified by two further measurements on the 
exchange process. From the previous two models, we can see that differ- 



factors over total money corresponding to their respec- 
tive holders after the steady distribution of money among 
agents had been observed. The value is 0.86 which certifies 
again that there is more money in the hands of holders 
with higher saving factors. The average value of saving 
factors over the money transferred was also computed, its 
value is about 0.52. This fact says that the money held 
by agents with higher saving factor has lower probability 
to take part in trade. If all money has equal probability, 
combining with the fact that the agents with higher sav- 
ing factors hold more money, it can be deduced that the 
value of this kind of average saving factor should be about 
0.86 all the time. Thus, we can conclude that the higher 
the saving factor of a unit of money's holder, the smaller 
probability for it to be transferred. 



5 Model with Preferential Spending 



Firstly, we measured the correlation coefficient of agents' ent saving patterns lead to different holding time distri- 

saving factors and the amount of money in their hands, butions. Especially, when the agents' saving factors are 

As shown in Figure 5, the correlation coefficient increases diverse, the probability of money to take part in trade 

sharply at the beginning of simulation, starts to decrease differs. Nevertheless, the probabilities of the money held 

slowly after about 2000th round. The reason of the re- by the same one agent are equal to each other. This is 

duction is that the correlation coefficient can not pick up an implicit assumption in all the simulations which means 
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the money is homogenous to any agent. However, it is not The power distribution arises from the diversity of the 
the case in real life. As the medium of exchange, money probability of money to participate in exchanges. At be- 
changes hand to hand. In this circulation process, money ginning of our simulation, no money has ever taken part 
abrades unavoidably. And when agents make exchange, in the trade, thus the probabilities are equal for all money 
the payers might spend their money with preference ac- according to Equation Q. After some of money are ex- 
cording to the degree of abrasion. As a result, the money changed randomly, they have higher probabilities and the 
is not homogenous for agents. To overcome this unrealistic others have relative lower ones. As the times of exchange 
feature, we proposed a new model which is quite similar increase, this slight diversity of money in the probability 
to the model with uniform saving factor. The only alter- will be enlarged till a stable distribution is formed. To see 
ation is that the probability of money chosen to change this process from another point of view, the longer for one 
hand is not equal even if the money is held by the same unit of money to wait, the lower probability for it to be 
agent, in other words, that the agents spend money with spent. In this way, comparing with the case without prefer- 
preference. ence, some money's holding times get shorter, while some 

In each round, two agents, i and j, are chosen ran- get much longer. Thus the power distribution appears, 

domly to participate in the trade. The amount of money We studied the holding time distribution at different 

transferred is determined by Equation (jSJ). If agent i is times, and found the power distribution is robust. For in- 

the payer, the probability of money k among m, to be stance, the holding time distribution for s = still has 

transferred is given by: the power form even after t = 500000. Contrarily to this, 

I _l_ i it is just after t — 1000 that one can clearly observe the 

p( k ) = —, ; ( 7 ) 



\^(l n + 1) steady distribution. 



71=1 

where l n is the times that money n has participated in the 
trade since the beginning of simulation. Here, we express 
the probability with the sum of exchange times and 1 in- In this paper, the effects of saving and spending patterns 



6 Conclusions 



stead of exchange times itself in case that denominator be on the distribution of money over holding time are ex- 
zero at the beginning of simulations. amined by computer simulations. All the simulations are 
The probability distributions of holding time for sev- performed basing on the ideal gas-like models. We consider 
eral different saving factors are recorded after money dis- two kinds of assumptions on saving pattern, one is that all 
tributions reach stationary state which are shown in Fig- agents have uniform saving factor, the other one is that 
ure 6. All distributions obey power law, and the only dif- the saving rates are set randomly distributed among the 
ference is the exponent. agents. In the model with uniform saving factor, the distri- 
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bution of money over the holding time follows an exponen- sium, edited by H. Takayasu (Springer, Tpkyo, 2003), 

tial law, while in the model with diverse saving factor the |Arxiv:cond -mat/0302147 

distribution changes to a power type. We further propose 8. A. Chatterjee, B. K. Chakrabarti and S. S. Manna, arXiv: 

a new trading model where the agents spend money with |cond-mat/0301289"t2003) . 

preference and also get power distribution. The simulation 9 ' A ' Chakraborti, Int. J. Mod. Phys. C. 13 1315-1321 (2002). 

• i , • ii 10. A. Dragulescu and V. M. Yakovenko, in Modelinq of Com- 
results indicate that the final distribution is determined by 

the character of the probability that money is chosen to 

participate in the trade. 



plex Systems: Seventh Granada Lectures, AIP Conference 
Proceedings 661 (American Institute of Physics, New York, 
2003), p. 180-183. 
11. A. Dragulescu, and V. M. Yakovenko, Physica A 299, 213 
(2001). 
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Fig. 1. Schematic presentation of the sampling method of 
holding time adopted: (a) in this paper; (b) in Reference |16|. to 
denotes the sampling time point, the light horizontal solid lines 
represent the evolution history of money, the vertical short bars 
symbolize moments for corresponding money to be transferred 
and then the dark segments correspond to the holding times 
to be recorded. 
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Fig. 2. The stationary holding time distribution obtained from 
the basic ideal gas-like model simulations versus holding time. 
The fitting in the inset indicates the distribution follows the 
exponential law: Ph{t) = ^ea;p(— r/T). 
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Fig. 3. The stationary distributions of holding time for several 
saving factors from to 0.9 derived from the simulations of the 
model with uniform saving factor in the semi-logarithmic scale. 
Note that in the figure the probabilities have been scaled by 
the maximum probability respectively. 
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Fig. 4. The stationary distribution of holding time derived 
from the simulations of the model with diverse saving factor 
in double logarithmic scale. The solid line is numerically fitted 
line in the form of Ph(t) oc t~ 1,14 . 
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Fig. 5. Correlative coefficient between the amount of money 
held by agents and their saving factors versus time. At time 
t — f54f , the coefficient reaches its maximum 0.773. 
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Fig. 6. The stationary distributions of holding time for several 
uniform saving factors from to 0.9 derived from the simu- 
lations of the model with preferential spending in the double 
logarithmic scale. Note that in the figure the probabilities have 
been scaled by the maximum probability respectively. 
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As well known, the exploration of the distribution of 
money can be traced back at least a century to the work 
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nomical and financial problems 1,2,3 . Recently, Some 
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might be some physical and mathematical rules govern- intervals for money to be held. This kind of time interval 
ing the distribution of income or wealth in the world and was ever called by Wicksell the "average period of idle- 
attempt to discover them. A series of models have been ness" or "interval of rest" of money |15| . In our previous 
developed for the equilibrium money distribution based on work , we called it "holding time" of money and found 
the analogy between market economics and kinetic theory that after the economy has achieved an equilibrium state, 
of gases 4,5 6,812 . Identifying exchange between any there is not only a distribution of money among agents, 
two agents in a closed economy where the money is con- but also a steady distribution over the holding time. We 
served with the two-body elastic collision in an ideal gas, also found that monetary velocity, an important macroe- 
these models show no matter how uniformly and forcefully conomic variable, which is associated with Irving Fisher 
one distributes money among agents initially, the succes- J7| , could be expressed as the expectation of the recipro- 
sive tradings eventually lead to a steady distribution of cal of holding time, 
money. And the shape of money distribution is determined 
by the trading rule for choosing an amount of money to 
transfer. Allowing agents to hold back some of their money 
when they are chosen to trade, B.K. Chakraborti et al. 

In a basic ideal gas-like model, the distribution of money 

introduced the saving behavior into the model by adding 

over the holding time follows an exponential law, where 

a saving factor s in the trading rule UJ. The simulation 

saving behavior is not taken into account. The purpose 

results clearly indicate a robust Gibbs-like distribution 
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a 'critical' Pareto distribution of money is found when 

change to power distribution. Then we introduce prefer- 

saving factor is set diversely among agents [E1EIIH]- 

ential spending behavior into the model in section [S] and 
In practice, money is held as a store of value, what again obtain power distribution. Comparing these results, 
is more, it plays an essential role for being a medium of we can conclude that the formation of holding time distri- 
exchange. Money is transferred consecutively from hand bution is associated with the character of the probability 
to hand in the exchange process, in which there exist time of money to be exchanged. 
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2 An Ideal Gas-Like Market Model and 
Holding Time Distribution 

We begin with the basic ideal gas-like model which was 
introduced firstly by A. Dragulescu and V.M. Yakovenko 
A close economy is considered in the model where 
the amount of money M is conserved and the number 
of agents N is fixed. The money is possessed by agents 
individually and agents can exchange money with each 
other. Since the scale and initial distribution of money 
have no effect on the final results, most of our simulations 
were carried out with N = 250 and M = 25000 and the 
amount of money held by each agent was set to be M/N 
at the beginning. The trade in the economy is modelled to 
take place round by round. In each round, two agents, i 
and j for example, are chosen randomly to get engaged in 
a trade among which agent i is "receiver" and the other 
one j is "payer" . The amount of money that changes hand 
Am is determined by trading rule which ensures that the 
amount of any agent's money is non-negative and the total 
money is conserved. A trading rule commonly used can be 
expressed as Am = s(rm + rrij)/2, where e is a random 
number from zero to unity. As for which units of money 
are chosen to be transferred, all in the payer's hand is 
equally probable. 

In the ideal gas-like model, money is held by agents 
and transferred frequently. In this process, if an agent re- 
ceives money from other agents, he will hold it in hand 
till paying it out to some other agents. The time interval 
between the receiving and paying out is named as holding 
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time |lfij . The holding times of a certain unit of money at 
different moments or those of different units of money at a 
given moment are not the same. We introduce the proba- 
bility distribution function of holding time Ph(r), which is 
defined so that the amount of money whose holding time 
lies between r and r + dr is equal to MPh{r)dr. So, we 
can get the normalization condition and the expression of 
the expectation of holding time as follows: 

POO 

/ P h {r)dT = l (1) 
Jo 

and 

/>oo 

T= / rP h (r)dT. (2) 
Jo 

In the simulations, suppose it is at round to that we 
start to reord, and so, holding time is recorded as the 
difference between the moments when the money takes 
part in trade after to for the first two times. The recording 
mode is illustrated in Figure 1(a). Please note this mode 
is different from what we have adopted in Reference |16|. 
which is shown in Figure 1 (b) . The measurement results of 
the two modes seem quite different, however, they reflect 
the same process in different ways. We adopt herein the 
mode (a) solely to facilitate the exposition. The typical 
distribution of holding time is shown in Figure 2. It can 
been seen from the inset of Figure 2 that the distribution 
of holding time follows an exponential law: 

iMr) = ie-*. (3) 

This result indicates that the transferring process of money 
is a Poisson process with intensity of ^. 

To get the distribution of holding time without system- 
atic factor disturbing, we performed the simulations about 
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100 times with different random seeds and data were not and the number of agents is fixed. Any agent's money 
collected until the probability distribution of money got is non-negative or no debt is allowed. The agents are in- 
stationary. And for convenience we stopped data collecting distinctive at the beginning of simulations: same initial 
after majority of money(> 99.9%) had been recorded. In amount of money and same saving factor s. In each round, 
all the following simulations, holding times are measured an arbitrary pair of agents are chosen to make exchange 
in this way after the distributions of money get stationary with each other. For example, at t-th round, agent i and j 
of course. take part in trading, so that at t + 1-th round their money 

rrii(t) and mj(t) change to 

3 Model with Uniform Saving Factor , , 

m;(f + 1) = m t {t) + Am; m,j(t + 1) = m 3 -(f) - Am; (4) 

In reality, saving behavior is a natural action pattern for 

where 

any economic agent. In order to insure future consump- 
tion, people always keep a part of their money as saving. ^ m = — s )i( £ ~ + £m j w]i (5) 
The ratio of the saving to total amount of money held by 

stitution, it is obvious that the trading rule satisfies the 



and e is a random fraction. After a straight-forward sub 

an agent is called "marginal propensity to save" by B. K 



Chakrabarti's group. The term of marginal propensity to . . ,. . .... . . 

conservation and non- negativity condition, and each agent 

save has totally different meaning in economics, which is . r , . . . . 

saves fraction s of his money before trade. 

defined as the partial derivative of saving function with re- _,. . . . . . . _,. „ „ 

Ihe simulation results are shown m figure 6, for some 

spect to income To avoid confusion, we rename it as . „ T , , , .. . . .. , ., 

^* values of s. It can be seen that the probability distribu- 

"saving factor". Referring to the saving factor, two cases . (.•■•,,. • r n n i 

tions of holding time for all saving factors decay cxponcn- 

have been considered by Chakrabarti's group, one is that . . , . . . . , . 

tially. And the lower the saving factor is, the steeper the 

all agents have a uniform saving factor, the other one is .... . .... , . , . , r 

distribution curve. I hese results indicate this kind of sav- 

that saving factors are randomly distributed among agents . , . . , . _ . r , . 

mg behavior does not change the Poisson nature of the 

IBHSI. As mentioned above, they found the equilibrium dis- , . . . ., . , 

exchanging process, but its intensity. 

tributions of money among agents had remarkable differ- 
ent characters under such two assumptions. Along this 

4 Model with Diverse Saving Factor 

line, in this section and the next one we shall examine 

the impacts of the saving behavior on the distribution of In realistic economy, how much an agent saves depends 
holding time for the two cases respectively. on the economic situations he or she faces, and the saving 

All the assumptions of the above ideal gas-like model factor of course varies from agent to agent due to their 
do work in this model. The amount of money is conserved different conditions. To get closer to reality, this model 
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inherits all the assumptions and evolution mechanism of non-linear associations. We also found the correlation co- 

the previous model except that of uniform saving factor, efficient falls to and keeps at about 0.32 after 500000th 

Each agent's saving factor is initialized at the beginning round. Although the value of the correlation coefficient 

of simulations which distributes randomly and uniformly is not high enough, it still implies that the agents with 

within an interval to 1, and is fixed in the simulations, higher saving factors hold more money. 

Correspondingly, the trading rule Equation © changes to Secondly, we computed the average value of saving 



Am = (1 - s t )(e - l)m t {t) + (1 - s^em^t)- (6) 



where Si, Sj are the saving factors of agent i and j respec- 
tively. 

To our surprise, once the diverse saving factor is intro- 
duced into the model, as shown in Figure 4, the holding 
time distribution changes to obey a power law instead of 
an exponential law. This result indicates that the transfer- 
ring process of money in this model is not a Poisson pro- 
cess any more. The Poisson nature of the process is broken 
due to the loss of homogeneity of the money transferring. 
In the previous model, for any saving factor, the proba- 
bility of each unit of money to participate in exchanges 
at any round is equal because the saving factor is set to 
be uniform for all agents. On the contrary, in this model 
the transferring probability of money is not the same any 
more due to the diversity of the saving factors. This con- 
clusion was verified by two further measurements on the 
exchange process. From the previous two models, we can see that differ- 



factors over total money corresponding to their respec- 
tive holders after the steady distribution of money among 
agents had been observed. The value is 0.86 which certifies 
again that there is more money in the hands of holders 
with higher saving factors. The average value of saving 
factors over the money transferred was also computed, its 
value is about 0.52. This fact says that the money held 
by agents with higher saving factor has lower probability 
to take part in trade. If all money has equal probability, 
combining with the fact that the agents with higher sav- 
ing factors hold more money, it can be deduced that the 
value of this kind of average saving factor should be about 
0.86 all the time. Thus, we can conclude that the higher 
the saving factor of a unit of money's holder, the smaller 
probability for it to be transferred. 



5 Model with Preferential Spending 



Firstly, we measured the correlation coefficient of agents' ent saving patterns lead to different holding time distri- 

saving factors and the amount of money in their hands, butions. Especially, when the agents' saving factors are 

As shown in Figure 5, the correlation coefficient increases diverse, the probability of money to take part in trade 

sharply at the beginning of simulation, starts to decrease differs. Nevertheless, the probabilities of the money held 

slowly after about 2000th round. The reason of the re- by the same one agent are equal to each other. This is 

duction is that the correlation coefficient can not pick up an implicit assumption in all the simulations which means 
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the money is homogenous to any agent. However, it is not The power distribution arises from the diversity of the 
the case in real life. As the medium of exchange, money probability of money to participate in exchanges. At be- 
changes hand to hand. In this circulation process, money ginning of our simulation, no money has ever taken part 
abrades unavoidably. And when agents make exchange, in the trade, thus the probabilities are equal for all money 
the payers might spend their money with preference ac- according to Equation Q. After some of money are ex- 
cording to the degree of abrasion. As a result, the money changed randomly, they have higher probabilities and the 
is not homogenous for agents. To overcome this unrealistic others have relative lower ones. As the times of exchange 
feature, we proposed a new model which is quite similar increase, this slight diversity of money in the probability 
to the model with uniform saving factor. The only alter- will be enlarged till a stable distribution is formed. To see 
ation is that the probability of money chosen to change this process from another point of view, the longer for one 
hand is not equal even if the money is held by the same unit of money to wait, the lower probability for it to be 
agent, in other words, that the agents spend money with spent. In this way, comparing with the case without prefer- 
preference. ence, some money's holding times get shorter, while some 

In each round, two agents, i and j, are chosen ran- get much longer. Thus the power distribution appears, 

domly to participate in the trade. The amount of money We studied the holding time distribution at different 

transferred is determined by Equation (jSJ). If agent i is times, and found the power distribution is robust. For in- 

the payer, the probability of money k among m, to be stance, the holding time distribution for s = still has 

transferred is given by: the power form even after t = 500000. Contrarily to this, 

I _l_ i it is just after t — 1000 that one can clearly observe the 

p( k ) = —, ; ( 7 ) 



\^(l n + 1) steady distribution. 



71=1 

where l n is the times that money n has participated in the 
trade since the beginning of simulation. Here, we express 
the probability with the sum of exchange times and 1 in- In this paper, the effects of saving and spending patterns 



6 Conclusions 



stead of exchange times itself in case that denominator be on the distribution of money over holding time are ex- 
zero at the beginning of simulations. amined by computer simulations. All the simulations are 
The probability distributions of holding time for sev- performed basing on the ideal gas-like models. We consider 
eral different saving factors are recorded after money dis- two kinds of assumptions on saving pattern, one is that all 
tributions reach stationary state which are shown in Fig- agents have uniform saving factor, the other one is that 
ure 6. All distributions obey power law, and the only dif- the saving rates are set randomly distributed among the 
ference is the exponent. agents. In the model with uniform saving factor, the distri- 
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bution of money over the holding time follows an exponen- sium, edited by H. Takayasu (Springer, Tpkyo, 2003), 

tial law, while in the model with diverse saving factor the cond-mat/0302147 

distribution changes to a power type. We further propose 8. A. Chatterjee, B. K. Chakrabarti and S. S. Manna, arXiv: 

a new trading model where the agents spend money with |cond-mat/0301289 

preference and also get power distribution. The simulation 9 ' A ' Chakraborti, Int. J. Mod. Phys. C. 13 1315-1321 (2002). 

• i , • ii 10. A. Dragulescu and V. M. Yakovenko, in Modelinq of Com- 
results indicate that the final distribution is determined by 

the character of the probability that money is chosen to 

participate in the trade. 



plex Systems: Seventh Granada Lectures, AIP Conference 
Proceedings 661 (American Institute of Physics, New York, 
2003), p. 180-183. 
11. A. Dragulescu, and V. M. Yakovenko, Physica A 299, 213 
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Fig. 1. Schematic presentation of the sampling method of 
holding time adopted: (a) in this paper; (b) in Reference |16|. to 
denotes the sampling time point, the light horizontal solid lines 
represent the evolution history of money, the vertical short bars 
symbolize moments for corresponding money to be transferred 
and then the dark segments correspond to the holding times 
to be recorded. 
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Fig. 2. The stationary holding time distribution obtained from 
the basic ideal gas-like model simulations versus holding time. 
The fitting in the inset indicates the distribution follows the 
exponential law: Ph{t) = ^ea;p(— r/T). 
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Fig. 3. The stationary distributions of holding time for several 
saving factors from to 0.9 derived from the simulations of the 
model with uniform saving factor in the semi-logarithmic scale. 
Note that in the figure the probabilities have been scaled by 
the maximum probability respectively. 
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Fig. 4. The stationary distribution of holding time derived 
from the simulations of the model with diverse saving factor 
in double logarithmic scale. The solid line is numerically fitted 
line in the form of Ph(t) oc t~ 1,14 . 
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Fig. 5. Correlative coefficient between the amount of money 
held by agents and their saving factors versus time. At time 
t — f54f , the coefficient reaches its maximum 0.773. 
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Fig. 6. The stationary distributions of holding time for several 
uniform saving factors from to 0.9 derived from the simu- 
lations of the model with preferential spending in the double 
logarithmic scale. Note that in the figure the probabilities have 
been scaled by the maximum probability respectively. 



